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ABSTRACT 

We report results for the alignments of galaxies in the EAGLE and cosmo-OWLS 
hydrodynamical cosmological simulations as a function of galaxy separation (—1 ^ 

logio(r/[/i _1 Mpc]) ^ 2) and halo mass (10.7 ^ logi 0 (M 2 oo/[ h~ x M 0 ]) < 15). We focus 
on two classes of alignments: the orientations of galaxies with respect to either the directions 
to, or the orientations of, surrounding galaxies. We find that the strength of the alignment is 
a strongly decreasing function of the distance between galaxies. For galaxies hosted by the 
most massive haloes in our simulations the alignment can remain significant up to ^ 100 Mpc. 
Galaxies hosted by more massive haloes show stronger alignment. At a fixed halo mass, more 
aspherical or prolate galaxies exhibit stronger alignments. The spatial distribution of satellites 
is anisotropic and significantly aligned with the major axis of the main host halo. The major 
axes of satellite galaxies, when all stars are considered, are preferentially aligned towards the 
centre of the main host halo. The predicted projected direction-orientation alignment, e g+ (r p ), 
is in broad agreement with recent observations. We find that the orientation-orientation align¬ 
ment is weaker than the orientation-direction alignment on all scales. Overall, the strength of 
galaxy alignments depends strongly on the subset of stars that are used to measure the ori¬ 
entations of galaxies and it is always weaker than the alignment of dark matter haloes. Thus, 
alignment models that use halo orientation as a direct proxy for galaxy orientation overesti¬ 
mate the impact of intrinsic galaxy alignments. 

Key words: cosmology: large-scale structure of the Universe, cosmology: theory, galaxies: 
haloes, galaxies: formation 


1 INTRODUCTION 

Tidal gravitational fields generated by the formation and evo¬ 
lution of large-scale structures tend to align galaxies due to 
correlations of tidal tor ques in random gaussian fields (e.g. 
iHeavens & PeacocSl 19881) . Analytic theories have been developed 
to describ e these l ar ge-scale alignments (li near alignment theory; 
ICatelan. Kamionkowski & BlandforcfeOQll) . but these are only ap¬ 
plicable to low matter density contrasts (the linear regime of struc¬ 
ture formation) and do not account for drastic events such as merg¬ 
ers of structures, which may erase initial correlations. 
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To overcome these limitations, galaxy alignments 
have been studied_via N- body s imulations (see e.g. 


West. Villumsen & Dekel 1991; |Tormenl 19971: 

Croft & Metzler 

2000 

;l 

Heavens. Refregier & Hevmansl;2000l; find 

120021 iLee et al. 

2008 

Bett 2012j). The most common ansatz in 

such studies is 


that galaxies are perfectly aligned with their dark matter haloes 
and that one can therefore translate the alignments of haloes 
directly into those of the galaxies that they host. However, the 
observed light from galaxies is emitted by the baryonic component 
of haloes and hydro-dynamical simulations of galaxy formation 
have revealed a misalignment between the baryonic and dark 
matter components of hal oes (iDeason et al.l l201ll : iTenneti et al.1 
120141 : [Velliscig et ai1l2015l) . On spatial scales characteristic of a 
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galaxy, baryon processes (radiative cooling, supernova explosions 
and AGN feedback) play an important role in shaping the spatial 
distribution of the stars that constitute a galaxy. Specifically, the 
ratio between cooling and heating determines the way baryons 
lose angular momentum and consequently the way they settle 
inside their dark matter haloes. Furthermore, feedback from star 
formation and AGN can heat and dis place large quantities of gas 
and inhibit star formation (Springe!. Di Matteo & Hernouistl 2005 


lgel, U 1 Matteo & hlernquist ZUID; 

I 2005: |Pi Matteo et al J 2008; 


Di Matteo, Springel & Hernquis 


Booth & Schavel 120091 : McCarthy et all l201oh . These processes, 
which determine when and where stars form, may influence 
the observed morphology of galaxies and in turn their observed 
orientations. In hydrodynamical simulations of galaxy formation 
in a cosmological volume, such processes are modelled simul¬ 
taneously, leading to a potentially more realistic realization of 
galaxy alignments. The study of such models can unveil patterns 
that encode important information concerning both the initial 
conditions that gave rise to the large-scale structure, and the 
evolution of highly non-linear structures like groups and clusters 
of galaxies. 

Beyond their relevance to galaxy formation theory, galaxy 
alignments are a potential contaminant of weak gravitational lens- 
ing measurements. Although this contamination is relatively mild, 
it is a significant concern for large-area cosmic shear sur veys 
(Ijoachimi et alJl20l5 iKirk et~aPl2015l : iKiessling et alJ 120151 and 
references therein). Requirements for the precision and accuracy of 
such surveys are very challenging, as their main goal is to constrain 
the dark energy equation-of-state parameters at the sub-percent 
level. Weak lensing surveys are used to measure the effect of the 
bending of light paths of photons emitted from distant galaxies due 
to intervening matter density contrasts along the line of sight. The 
distortion and magnification of galaxy images is so weak that it can 
only be characterized by correlating the shapes and orientations of 
large numbers of background galaxies. In a pure weak gravitational 
lensing setting, the observed ellipticity of a galaxy, e, is the sum of 
the intrinsic shape of the galaxy, e s , and the shear distortion that 
the light of the galaxy experiences due to gravitational lensing, 7 , 


e + 7 • 


( 1 ) 


If galaxies are randomly oriented, the average ellipticity of a sample 
of galaxies, (e s ), vanishes. Therefore, any detection of a nonzero 
(e) is interpreted as a measurement of gravitational shear 7 . How¬ 
ever, in the limit of a very weak lensing signal, the distortion in¬ 
duced via gravitational forces (giving rise to an intrinsic alignment) 
can be a non-negligible fraction of the distortion due to the pure 
gravitational lensing e ffe ct (often termed appar ent alignment, see 
ICrittenden et al.l200ll and lCrittenden et al.ll2002l for a statistical de¬ 
scription of this effect). 

Cosmic shear measurements are obtained in the form of pro¬ 
jected 2-point correlation functions (or their equivalent angular 
power spectra) between shapes of galaxies. Following Eq.Q] 

( ee ) = ( 77 ) + ( 7 e s ) + © 7 ) + (eV) , (2) 

= GG + GI + IG + II. (3) 


If we assume that galaxies are not intrinsically oriented towards one 
another, then the only correlations in the shape and orientation of 
observed galaxies is due to the gravitational lensing effect of the in¬ 
tervening mass distribution between the sources and the observer, 
( 77 ). In this case the only nonzero term is the GG (shear-shear) 
auto correlation. In the case of a non negligible intrinsic alignment 
of galaxies, the II term is also nonzero, i.e. part of the correla¬ 
tion between the shape and orientation of galaxies is intrinsic. If 


the same gravitational forces that shear the light emitted from a 
galaxy also tidally influence the intrinsic shape of other galaxies, 
then this will produce a nonzero cross correlation between shear 
and intrinsic shape (GI). The term IG is zero since a foreground 
galaxy cannot be lensed by the same structure that is tidally influ¬ 
encing a background galaxy, unless their respective position along 
the line of sight is confused due to large errors in the redshift mea¬ 
surements. 

In this paper we report results for the intrinsic alignment of 
galaxies in hydro-cosmological simulations. Specifically, we fo¬ 
cus on the orientation-direction and orientation-orientation galaxy 
alignments. To this aim, we define as galaxy orientation the major 
eigenvector of the inertia tensor of the distribution of stars in the 
subhalo. We then compute the mean values of the angle between 
the galaxy orientation and the separation vector of other galaxies, 
as a function of their distance. In the case of orientation-orientation 
alignment we compute the mean value of the angle between the ma¬ 
jor axes of the galaxy pairs, as a function of their distance. While 
the orientation-orientation alignment can be interpreted straightfor¬ 
wardly as the II term in Eq. m, the orientation-direction isr elated 
to the GI term in a less direct way (see lloachimi et al.1120111 for a 
derivation of the GI power spectrum from the ellipticity correlation 
function). 

In this paper we make use of four complementary simulations 
to explore the dependence of the orientation-direction alignment 
over four orders of magnitude in subhalo mass, and spanning phys¬ 
ical separations of hundreds of Mpc. The use of four simulations 
of different cosmological volumes offer both resolution and statis¬ 
tics, whilst also incorporating baryon physics. The EAGLE simu¬ 
lations used in this work have been calibrated to reproduce the ob¬ 
served present-day galaxy stellar m ass functio n and the observed 
size-mass relation of disc galaxies (|Schave et all 20151). wh ereas 
the cosmo-OWLS jLe Brun et al.l20l4lMcCarthv et alJ2014h sim¬ 
ulations reproduce key (X-ray and optical) observed properties of 
galaxy groups and clusters, in addition to the observed galaxy mass 
function for haloe s more massive than log(M/[/i -1 M©]) = 13. 
In IVelliscig et al.l d2015h we used the same set of simulations to 
study the shape and relative alignment of the distributions of stars, 
dark matter, and hot gas within their own host haloes. One of the 
conclusions was that although galaxies align relatively well with 
the local distribution of the total (mostly dark) matter, they exhibit 
much larger misalignments with respect to the orientiation of their 
complete host haloes. 

Afte r t he sub mission of this manuscript, a paper by 
IChisari etal ]& appeared on the arXiv. They study the align¬ 
ment of galaxies at z = 0.5 in the cosmolo gical hydrodynami¬ 
cal simulation HORIZON-AGN (iDubois et al. 2014 ) r un wi th the 
adaptive-mesh-refinement code RAMSES (iTevssierl 12002 ). The 
HORIZON-AGN simulation is run in a (100 h~ x Mpc) 3 volume 
with a dark matter particle mass resolution of ma m = 8 x 
10 7 M©. They focus on a galaxy stellar mass range of 9 < 
logio(M sta r/[M©]) < 12.36 and separations up to 25 h~ x Mpc. 
Their analysis differs in various technical, as well as conceptual, 
aspects from the study presented here. However, they also report 
that the strength of galaxy alignments depends strongly on the sub¬ 
set of stars that are used to measure the orientations of galaxies, as 
found in our investigation. 

Throughout the paper, we assume a flat ACDM cosmology 
with massless neutrinos. Such a cosmological model is character¬ 
ized by five parameters: {Q m , Ob, cr 8 , n s , h}. The EAGLE and 
cosmo-OWLS simulations were run with two slightly different sets 
of values for these parameters. Specifically, EAGLE was run us- 
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Table 1 . List of the simulations used and their relevant properties. Description of the columns: (1) descriptive simulation name; (2) comoving size of 
the simulation box; (3) total number of particles; (4) cosmological parameters; (5) initial mass of baryonic particles; (6) mass of dark matter particles; 
(7) maximum proper softening length; (8) simulation name tag. 


Simulation 

(1) 

L 

(2) 

-^particle 

(3) 

Cosmology 

(4) 

m h 

[^M 0 ] 

(5) 

TO dm 

[fc _1 M©] 

(6) 

Cprop 

[ h~ x kpc] 
(7) 

tag 

(8) 

EAGLE Recal 

25 [Mpc] 

2 x 752 3 

PLANCK 

1.5 x 10 5 

8.2 x 10 5 

0.2 

EA L025 

EAGLE Ref 

100 [Mpc] 

2 x 1504 3 

PLANCK 

1.2 x 10 6 

6.6 x 10 6 

0.5 

EAL100 

cosmo-OWLS AGN 8.0 

200 [h- 1 Mpc] 

2 x 1024 3 

WMAP7 

8.7 x 10 7 

4.1 x 10 8 

2.0 

CO L200 

cosmo-OWLS AGN 8.0 

400 [h- 1 Mpc] 

2 x 1024 3 

WMAP7 

7.5 x 10 8 

3.7 x 10 9 

4.0 

CO L400 


ing the set of cosmological values suggested by the Planck mission 
{Q m , O b , cr 8 , ns, h} = 10.307, 0.04825 0.8288, 0.9611, 0.6777} 
(Table 9: lPlanck Collaboration et all2014l) . whereas cosmo-OWLS 
was run using the cosm ological pa rameters suggested by the 7th- 
year data release jKomatsu et alii Il201 lh of the WMAP mission 
{O m , Q b , cr 8 , n s , h} = {0.272, 0.0455, 0.728, 0.81, 0.967, 0.704}. 

This paper is organized as follows. In Section |3 we summa¬ 
rize the properties of the simulations employed in this study (§ 12.11) 
and we introduce the technical definitions used throughout the pa¬ 
per (§ !2.2l and § 12.31) . In Section [3] we report the dependence of the 
orientation-direction alignment of galaxies on subhalo mass (§ 13.11) . 
matter components (§ 13.21) . galaxy morphology (§ 13.31) and subhalo 
type (§ 13.41 ). In Section[4]we compare our results with observations 
of the orientation-direction alignment. In section[5]we report results 
for the orientation-orientation alignment of galaxies. We summa¬ 
rize our findings and conclude in Section[6] 


2 SIMULATIONS AND TECHNICAL DEFINITIONS 

2.1 Simulations 


In this work we employ two diffe rent sets of hydrody n amical cos¬ 
molo gical simulations, E AGLE (Schaveet al. 120151: Crain e t al] 
20151) and cosmo-OW LS dSchave et al]l20ld : lLe Brun et al.ll2014l : 


McCarthy et alJ 120141) . Specifically, from the EAGLE project 


we make use of the simulations run in domains of boxsize 
L — 25 and 100 comoving Mpc in order to study with suf¬ 
ficient resolution central and satellite galaxies hosted by sub¬ 
haloes with mass from logio(M su b/[ h~ x M©]) = 10.7 up to 
logio(M su b/[/i -1 M©]) = 12.6, whereas from cosmo-OWLS we 
select the simulations run in domains of boxsize L = 200 and 
400 comoving h~ x Mpc which enable us to extend our analysis to 
logio(M su b/[/i -1 M©]) = 15. For each simulation the minimum 
value of subhalo mass is chosen to be the subhalo mass above which 
all haloes have at least 300 stellar particles. Using 300 parti cles 
ensures a reliable estimation of the subhalo shape dVelliscig et alJ 
20151) . Table |T] lists relevant specifics of these simulations. A rel¬ 
evant feature of our composite sample of haloes, taken from four 
different simulations, is that it reproduces the stellar mass halo 
mass relatio n inferred from abundance matching techniques studies 
dSchave et al.l2015h . which ensures that galaxies in our simulations 
reside in subhaloes of the right mass. 

EAGLE and cosmo-OWLS were both run using modified 
versions of the iV-Body Tree-PM smooth ed par ticle hydrody¬ 
namics (SPH) code GADGET 3 JSpringell l2005h . The simula¬ 
tions employed i n this work make u se of e lement-by-element ra¬ 
diative cooling (Wiersma. Schave & Smitbll2009h . star formation 
dSchave & Dalla Vecchial2008b . stellar mass losses (I Wiersma et alJ 


12009), stellar feedback toalla Vecchia & Schavel I2008L 2012), 
Black Hole (BH) growth through gas accretion and mergers 
teooth & Schavel UOO^; iRosas-Guevara et al]l2013l). and thermal 
AGN feedback ( Booth & Schavell2009l : ISchave et al.ll2015l) . 

The subgrid physics used in cos mo-O WLS is identical to 
that used in the OWLS run ”AGN” (ISchave et al.112010 ). EA¬ 
GLE includes a series of developments with respect to cosmo- 
OWLS in the subgrid physics, namely the use of thermal 
flDalla Vecchia & Schavel 12012) . instead of kinetic, energy feed¬ 
back from star for mation, BH accretion that d epends on the gas an¬ 
gular momentum (IRosas-Guevara et aOl2013l) and a metallicity de¬ 
pendent star formation law. More information regarding the techni¬ 
cal implementation of EAGLE’s hydr o-dynamical a s pects, as well 
as the subgrid physics, can be found in lSchave et al.1 d2015h . 


2.2 Halo and subhalo definition 

Haloes are identified by first applying the Friends-of-Friends (FoF) 
algorit hm to the dark matter particles, with linking length 0.2 
toavis et al.ll 19851) . Baryonic particles are associated to their clos¬ 
est dark matter particle and they inherit their group classification. 
Subhaloes are identified as groups of particles in local minima of 
the gravitational potential. The gravitational potential is calculated 
for the different particle types separately and then added in order 
to avoid biases due to different particle masses. Local minima are 
identified by locating saddle points in the gravitational potential. 
All particles bound to a given local minimum constitute a subhalo. 
The most massive subhalo in a given halo is the central subhalo, 
whereas the others are satellite subhaloes. Minima of the gravita¬ 
tional potential are used to identify the centers of subhaloes. The 
subhalo mass M sub is the sum of the masses of all the particles be¬ 
longing to the subhalo. For every subhalo we define the radius r b ^f 
within which half the mass in dark matter is found. Similarly, but 
using stellar particles, we define f (usually around one order of 
magnitude smaller than r b ^fX which represents a proxy for the typ¬ 
ical observable extent of a galaxy within a subhalo. The 7*200 is the 
radius of the sphere, centered on the central subhalo, that encom¬ 
passes a mean density that is 200 times the critical density of the 
Universe. The mass within 7*200 is the halo mass Mg p . The afore¬ 
mentioned quantities are computed using SUB FIND (iSpringel et al.1 
l200lUDolag etai]l2009l) . 

In Table |2] we summarize the z = 0 values of various quanti¬ 
ties of interest for the halo mass bins analysed here. 


2.3 Shape parameter definitions 

To describe the morphology and orientation of a subhalo we make 
use of the three-dimensional mass distribution tensor, also referred 
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Table 2. Values at z = 0 of various quantities of interest for our four subhalo mass bins. Description of the columns: (1) simulation tag; (2) 
subhalo mass range log 10 (M su b/( h _1 Mq)); (3) median value of the halo mass logio(M|J 1 0 t ) for centrals; (4) median value of the stellar mass 
(log 10 (M s tar/( h~ x Mq))); (5) standard deviation of the stellar mass distribution criogi 0 M st ar ’ (6) median value of halo virial radius t^qq for 
centrals; (7) median radius within which half of the mass in dark matter is enclosed; (8) median radius within which half of the mass in stars is 
enclosed; (9) number of haloes; (10) number of satellite haloes (11) color used throughout the paper for this particular mass bin and, with different 
shades, for the simulation from which the mass bin is drawn from. 


Simulation tag 

(1) 

mass bin 

* 

(2) 

* 

(3) 

Afstar 

* 

(4) 

^logioMstar 

* 

(5) 

-crit 
' 200 

** 

(6) 

^dm 
' half 
** 

(7) 

^ star 
' half 
** 

(8) 

Ahalo 

(9) 

77sat 

(10) 

Color 

(11) 

EA L025 

[10.70 - 11.30] 

10.87 

8.72 

0.46 

68.1 

28.0 

2.3 

234 

43 

black 

EAL100 

[11.30 - 12.60] 

11.59 

9.92 

0.45 

118.4 

50.7 

3.2 

4530 

745 

red 

CO L200 

[12.60 - 13.70] 

12.78 

10.88 

0.27 

295.6 

175.7 

31.1 

5745 

450 

green 

CO L400 

[13.70 - 15.00] 

13.82 

11.85 

0.22 

656.3 

416.4 

73.5 

3014 

94 

blue 


*logi 0 [M/(h 1 M 0 )] 
** [ h~ x kpc] 


to as the inertia tensor ( e.g. lCole & Lacevll 1996b . 


-Wpart 

Mij — ^ ^ TYlpXpiXpj , ( 4 ) 

P =1 

where iV part is the number of particles that belong to the structure 
of interest, x P i denotes the element i (with i,j = 1, 2, 3 for a 3D 
particle distribution) of the position vector of particle p, and m p is 
its mass. 

The eigenvalues of the inertia tensor are A^ (with i = 1, 2,3 
and Ai > A 2 > A 3 , for a 3D particle distribution as in our case). 
The moduli of the major, intermediate, and minor axes of the el¬ 
lipsoid that have the same mass distribution as the structure of in¬ 
terest, can be written in terms of these eigenvalues as a — y/Xi, 
b = \/A 2 , and c = VX 3 . Specific ratios of the moduli of the 
axes are used to define the sphericity, S = c/a, and tria xiali ty, 
T = (a 2 — b 2 )/(a 2 — c 2 ), parameters (see I Velliscig et al.|[2015h . 
The eigenvectors e*, associated with the eigenvalues Ai, define the 
orientation of the ellipsoid and are a proxy for the orientation of 
the structure itself. We interpret this ellipsoid as an approximation 
to the shape of the halo and the axis represented by the major eigen¬ 
vector as the orientation of the halo in a 3D space. 


3 ORIENTATION-DIRECTION ALIGNMENT 



ei ~ 

Orientation sample ( + ) 


Figure 1. Diagram of the angle <j> between the major eigenvector ei of a 
subhalo in the orientation sample (+), and the separation vector d pointing 
towards the direction of a subhalo in the position sample (g). Note that all 
quantities are defined in 3D space. 


In this section we present results concerning the alignment between 
the orientations of the stellar distributions in subhaloes, defined as 
the major eigenvector of the inertia tensor, ei, and the normalized 
separation vector, d , of a galaxy at distance r. Note that all quanti¬ 
ties are defined in a 3D space. We define 4> as: 

<j)(r) — arccos(|ei • d(r )|), (5) 

where ei is the major eigenvector of a galaxy in the orientation 
sample, and d is the separation vector pointing towards the position 
of a galaxy in the position sample (see Fig. |TJ. Note that, follow- 
ing Eq|5] 0 < <f> < 7T/2. The value of (cos(0)) is then computed 
as an average over pairs of galaxies from the orientation and po¬ 
sition samples. Values of (cos(0)) close to unity indicate that on 
average galaxies are preferentially oriented towards the direction 
of neighbouring subhaloes. We remind the reader that we use the 
term subhalo to refer to the ensemble of particles bound to a local 
minimum in the gravitational potential. Central galaxies are hosted 


by the most massive subhalo in a FoF group (see § 12.21 ). Through¬ 
out the text and in the figures we use (+) to refer to properties of 
galaxies in the orientation sample, whereas we use (g) for galaxies 
in the position sample. 

Observations typically measure the product of the cosine of 
the angle <j) and the ellipticity of the galaxy in the orientation sam¬ 
ple. We opt to begin our analysis by presenting results only for the 
angle cj) since it has a clearer interpretation that is independent on 
the shape determination of the galaxy. We present results for obser- 
vationally accessible proxies in Section |4] 

3.1 Dependence on subhalo mass and separation 

The left panel of Fig. |2] shows (cos(0)) for pairs of galaxies (both 
centrals and satellites) binned in subhalo mass and as a function of 
3D separation r. Subhaloes in the orientation sample (+) are cho- 
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Figure 2. Left : Mean value of the cosine of the angle f between the major eigenvector of the stellar distribution and the directions towards subhaloes with 
comparable masses as a function of 3D galaxy separation. Every mass bin is taken from a different simulation. The simulation identifiers used in the legends 
refer to column (8) of Table [I] The minimum subhalo mass in every bin ensures that only haloes with more than 300 stellar particles are selected. The curves 
are not shown for 3D separations larger than approximately 1/3 of the simulation volume. Right: Same as left panel but with physical distances rescaled by 
the f of the subhaloes. In both panels the error bars represent one sigma bootstrap errors. The horizontal dashed line indicates the expectation value for 
random orientations. The orientation-direction alignment decreases with distance and increases with mass. The mass dependence is greatly reduced when the 
distances are normalized by r^ f . 




Figure 3. As for the right panel of Fig. [2] but in this case the masses of the subhaloes in the orientation sample (+) are kept fixed whereas subhaloes in the 
position (g) sample are selected from mass bins above, below or equal to the mass bin of the orientation sample. Physical distances are rescaled by the of 
the subhaloes in the orientation sample. In the left panel the subhaloes are taken from the EAGLE LI00 simulation and in the right panel they are taken from 
the cosmo-OWLS L200 simulation. In both panels the error bars represent one sigma bootstrap errors. Thicker lines indicate more massive subhaloes for the 
position sample. The orientation-direction alignment is stronger for more massive subhaloes in the position subsample. 


sen to have the same mass limits as the subhaloes in the position 
sample (g). Values are shown for four different choices of subhalo 
masses, where every mass bin is taken from a different simula¬ 
tion (see legend). Errors are estimated via the bootstrap technique. 
Specifically, we use the 16th and the 84th percentiles of 100 re¬ 
alizations to estimate the lower and upper limits of the error bars. 
The cosine of the angle between the orientation of galaxies and the 


direction of neighbouring galaxies is a decreasing function of dis¬ 
tance and it increases with mass. For large separations the angle 
tends to the mean value for a randomly distributed galaxy orienta¬ 
tion, i.e. (cos ((/))) = 0.5. The physical scale at which this asymp¬ 
totic behaviour is reached increases with increasing subhalo mass. 

The right panel of Fig. |2] shows (cos(0)) as a function of 
the physical separation rescaled by the average size of subhaloes, 
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( r h"if}> in that mass bin. This rescaling removes most, but not all, 
of the offset between different halo mass bins. On average, subhalo 
pairs separated by more than lOOr^f show only weak alignment 
((cos(0)) ^ 0.52 at 100 r b ™ f ). 

Fig.[3]shows (cos(0)) as a function of the separation rescaled 
by the average size of the subhaloes in the orientation (+) sample. 
In this case the masses of the subhaloes for which we measure the 
orientation of the stellar distribution are kept fixed whereas haloes 
in the position (g) sample are selected from mass bins above, be¬ 
low or equal to the mass bin of the orientation sample. Results are 
shown for two of the four simulations: in the left panel for EAGLE 
LI00 and in the right panel for cosmo-OWLS L200. The line thick¬ 
ness is proportional to the subhalo mass of the position sample. The 
orientation of subhaloes of a given mass tends to be more aligned 
with the position of higher-mass subhaloes. 

We note that the two suites of simulations employed here, 
cosmo-OWLS and EAGLE, differ in resolution, volume, cosmol¬ 
ogy and subgrid physics. Testing how each of these differences im¬ 
pacts our mean results is beyond the scope of this study (we would 
need as many simulations as differences that we wish to test), there¬ 
fore we examine the overall convergence of the two simulations by 
selecting a subhalo mass bin 12.6 < logio(M sub /[/i 1 M®]) < 
13.1 that yields an orientation sample of galaxies that is numerous 
enough in the EAGLE LI00 simulation, as well as resolved in the 
cosmo-OWLS L200 simulation. We find that, in this specific case, 
the results are consistent within the bootstrapped errors, both for 
stars and stars within r^jf (not shown). 

Subhalo mass plays an important role in the strength of the 
orientation-direction alignment of subhaloes. The dependence on 
the subhalo mass weakens with distance but only becomes negligi¬ 
ble for separation ^>100 times the subhalo radius. 

3.2 Dependence on the choice of matter component 

In this section we report the orientation-direction alignment for the 
case in which the orientation of the subhalo is calculated using, re¬ 
spectively, dark matter, stars (as in the previous section) and stars 
within the half-mass radius • An alternative choice of a proxy 
for the typical extent of a galaxy would be to consider only stars 
within a fixed 3D aperture of 30Kpc that gives similar galaxy prop¬ 
erties a s the 2-D Pe trosi an apertures often used in observational 
studies dSchave et al.ll201 5). Note that the two definitions coincide 
for the subhalo mass bin 12.6 < logio(M su b/[/i _1 M 0 ]) < 13.7 
(CO L200). We note that 7*^5 varies among the four mass bins 
used in this work (see Table [2] column ( 8 )). 

Fig-Elshows the cosine of the angle cj) between the direction of 
nearby subhaloes and the orientation of the distribution of dark mat¬ 
ter, stars (as shown in Fig. [2]) and stars within the half-mass radius 
of subhaloes in the same mass bin. The left panel displays the re¬ 
sults for the subhalo mass bin 11.3 < logio(M sub /[(i 'Mej) < 

12.6 (from the EAGLE LI00 simulation), whereas the right panel 
refers to the subhalo mass bin 12.6 < logio(M su b/[ /i _1 M 0 ]) < 

13.7 (from the cosmo-OWLS L200 simulation). 

Irrespective of the subhalo mass and separation, the orienta¬ 
tion of the dark matter component shows the strongest alignment 
with the directions of nearby haloes, whereas the orientation of 
stars inside shows the weakest alignment. 

These results are suggestive of a scenario in which the align¬ 
ment between subhaloes and the surrounding density field is im¬ 
printed mostly on the dark matter distribution. Therefore, when 
the orientation of the subhalo is computed using all stars or the 
stars within f, the signal is weakened according to the internal 


misalignment angle between the specified component and the to¬ 
tal dark matter distribution. The trend shown by Fig. 0] therefore 
follows naturally from the results of IVelliscig et alJ ( 2015 ): stars 
within r^f exhibit a weaker alignment with the total dark matter 
distribution than all stars in the subhalo. 

The difference between the orientation-direction alignment 
obtained using the dark matter, all the stars or the stars within 
the typical extent of the galaxy, could account for the com¬ 
mon finding reported in the literature of galaxy alignment, that 
such alignments are systematically stronger in simulations than 
when measured in o bser vational data (se e the recent reviews of 
iKiessling et al.l 120151 and iKirk et all 120151 for a detailed compar¬ 
ison between observational and computational studies). Observa¬ 
tions are limited to the shape and orientation of the region of 
a galaxy above a limit surface brightness (often within surface 
brightness isophotes), whereas simulations need to rely on prox¬ 
ies for the ext ent of those regions (e.g. using baryonic overden¬ 
sity thresholds lHahn 1 Tevssier & Carolloll 2010 l : ICodis et aPl2015l : 
IWelker et~akl l20l4 iDubois et al.l l2014l) or to employ weighting 
schemes to the sample of star particles that constitute a gal axy (see 
e.g. use of the reduced inertia tensor in lTenneti et al .1120151) . 

3.3 Dependence on galaxy morphology 

Theory predicts that the alignment of early-type galaxies and 
late-type galaxies arises from differ ent p hysical processes (e.g. 
ICatelan. Kamionkowski & Blandfordll200ll) . It is of interest then 
to study the alignment as a function of galaxy morphologies. 

In this section we report the orientation-direction alignment 
of galaxies with different sphericities in order to explore the effect 
of the shape of galaxies on the orientation-direction alignment. We 
divide our sample of subhaloes according to the sphericity of their 
whole stellar distribution, defined as S = c/a where a and c are 
the squareroot of the major and minor eigenvalues of the inertia 
tensor respectively (see iO. We choose a threshold value for the 
sphericity of 0.5 that yields a similar numbers of galaxies in the two 
subsamples, as the median sphericity of the total sample is 0.55. 
This galaxy selection by sphericity represent a simple proxy for 
galaxy morphology. 

Fig. a shows the mean values of the cosine of the angle 0 
for galaxies of sphericity above and below the threshold, as well 
as for the total sample. The left panel displays the results for the 
subhalo mass bin 11.3 < logio(M su b/[ h~ x M 0 ]) < 12.6 (from 
the EAGLE LI00 simulation), whereas the right panel refers to the 
subhalo mass bin 12.6 < logio(M su b/[ h~ x M 0 ]) < 13.7 (from 
the cosmo-OWLS L200 simulation). 

More spherical galaxies (thinner lines) show a weaker 
orientation-direction alignment. The differences between the two 
shape selected samples of haloes are within the errors for scales 
larger than 1 h~ x Mpc, suggesting that the effect of shape is dom¬ 
inated by subhaloes of the same hosts. A similar trend (not shown) 
is found using triaxiality, see T2.31 as the indicator of galaxy 
shape. Prolate (T > 0.5) stellar distributions show the strongest 
orientation-direction alignment, whereas oblate (T < 0.5) ones 
show the weakest. The better alignment of prolate or aspherical 
galaxies is probably due to the fact that these galaxies align bet¬ 
ter with their underlying dark matter distributions (not shown), 
which in turn produces a stronger orientation-direction alignment 
(see Fig. s. 

We note that the orientation of a perfectly spherical distribu¬ 
tion (S = 1) of stars is ill defined. Although this can potentially af¬ 
fect our measurements, less than 2 % of galaxies in our sample have 
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r [Mpc/h] 



Figure 4. Mean value of the cosine of the angle (f> between the major eigenvectors of the distributions of stars (red curve in the left panel and green curve 
in right panel as in Fig. [2j, dark matter (gray curves), or stars within (purple curves) and the direction towards subhaloes with comparable masses as 
a function of 3D galaxy separation. The subhaloes used for the left panel are taken from the EAGLE LI00 (11. 3 < logio(M sub /[/i -1 M©]) < 12.6) 
simulation while in the right panel they are taken from the cosmo-OWLS L200 simulation (12.6 < logio (Mmb / [ h~ 1 Mq] ) < 13.7). Thicker lines indicate 
components with stronger alignment. In both panels the error bars represent one sigma bootstrap errors. The orientation of the dark matter component is most 
strongly aligned with the directions of nearby subhaloes, whereas the orientation of stars inside shows the weakest alignment. 




r [Mpc/h] r [Mpc/h] 


Figure 5. Mean value of the cosine of the angle (f> between the major eigenvector of the stellar distribution and the direction towards neighbouring subhaloes 
as a function of 3D galaxy separation, for galaxies in the orientation sample selected based on their shape. The selection is based on the sphericity of the whole 
stellar distribution defined as S' = c/a where a and c are the square root of the major and minor eigenvalues of the inertia tensor respectively. We choose a 
threshold value for the sphericity of 0.5. The subhaloes used for the left panel are taken from the EAGLE LI00 (11.3 < logio(M sub /[/i _1 Mq]) < 12.6) 
simulation while in the right panel they are taken from the cosmo-OWLS L200 simulation (12.6 < logio (M su b/ [ Mq]) < 13.7). Thicker lines 
indicate components with stronger alignment. In both panels the error bars represent one sigma bootstrap errors. More spherical galaxies show a weaker 
orientation-direction alignment. 


a sphericity higher than 0.8. We also note that more massive haloes, 
for which the orientation-direction align ment is strongest, tend to 
be less spherical and more triaxial fsee lVelliscig et al.l2015h . There¬ 
fore, selecting haloes by shape biases the sample towards system¬ 
atically different masses: however, the mass difference in the two 


shape-selected samples is about 4%, which is too small to explain 
the differences in alignment of haloes with different shapes. 

Observations indicate that ellipsoidal gal axies show s tronger 
intrinsic alignment than blue disk galaxies dHirata et al.l 120071 : 
ISingh. Mandelbaum & Morel 120151) . However, we caution the 
reader that there are still many complications to take into account 
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before one can compare the trends discussed above with these ob¬ 
servational results. First, one would need to select galaxies based 
on their colors, which requires stellar population synthesis mod¬ 
els. Second, the sphericity of the stellar component is a simplistic 
proxy for selecting disc galaxies. Selecting galaxies according to 
their morphology, in a similar way as done observationally, would 
require a stellar light decomposition in bulge and disc component. 


3.4 Alignment of satellite and central galaxies 

3.4.1 The increased probability of finding satellites along the 
major axis of the central galaxy 

In the previous sections we studied the orientation-direction 
alignment of galaxies irrespective of their classification as cen¬ 
trals or satellites. In this subsection we report the alignment 
between the orientations of central galaxies (g) and the di¬ 
rections of satellit43 galaxies (+) and, in turn, the probabil¬ 
ity of finding satellite galaxies distributed along the major axis 
of the central galaxy. This effect has been studied both theo- 


retically, making use of N-Body (e.g. 

Faltenbacher et al.1 120081: 

lAgustsson & Brainerd 20ld:lWan 

g et al. 

2014a) and hydrodynam- 

ical simulations jLibeskind et al. 

I 2007 I 

Deason et al. 

[2011), and 

observationallv dSales & Lambas 

2004; 

Brainerd 2005 

uYang et al. 

200^ Wang et al.ll2008:llNierenberg et alfep 121:11 Wang et alJ2014bl: 

Dong et al. 2014k Those studies report that the distribution of satel- 


lites around central galaxies is anisotropic, with an excess of satel¬ 
lites aligned with the major axis of the central galaxy. 

Fig.lH shows the average angle between the orientation of 
the stellar distribution of central subhaloes and the position of 
satellite galaxies. Values of (cosfi) that are significantly greater 
than 0.5 indicate that the positions of satellites are preferen¬ 
tially aligned with the major axis of the central galaxy. We use 
two different mass bins taken from two simulations: 11.3 < 
logioCMsub/l^M©]) < 12.6 from EAGLE L100 (left) and 
12.6 < logio(M su b/[/i -1 M©]) < 13.7 from cosmo-OWLS 
L200 (right). The line thickness is proportional to the subhalo mass 
of the position (g) sample. In both panels the physical separations 
between the pairs are normalized by the ( 7 * 200 ) °f the haloes host¬ 
ing the central galaxies. 

For separations up to 100 (r 200 )> the positions of satellite 
galaxies are significantly aligned with the orientation of central 
galaxies (not necessarily in the same host halo), with more massive 
satellites showing a stronger alignment. The same qualitative be¬ 
haviour is found for both mass bins, but the effect is stronger for the 
more massive central subhaloes. On scales larger than ~ (10 7 * 200 ) 
the alignment depends only weakly on the mass of the satellite sub¬ 
haloes. We speculate that the alignment of satellites with central 
galaxies of different host haloes is likely driven by the correlation 
between the orientation of the central galaxies and the surround¬ 
ing large-scale structure, which in turn influences the positions of 
satellite galaxies. 


3.4.2 The radial alignment of satellite galaxies with the direction 
of the host galaxy 

Here we investigate the radial alignment of the orientations of satel¬ 
lites (+) with the direction of the central galaxy (g), whereas in the 


previous section report the results for alignment between the ori¬ 
entations of the central galaxy and the direction of satellites. The 
orientation of satellite subhaloes is computed using all the stars 
bounded to the subhalo. Theoretical studies using N-body sim¬ 
ulatio ns ( Kuhlen. Diemand & Madaul [20071 : IPereira. Brvan & Gilll 
120081: Faltenbacher et al.l 120081) and hydrodynamic simulations 
llKnebe et alJl2010h found that on average the orientation of satel¬ 
lite galaxies is aligned with the direction of the centre of their host 
halo. 


Fig. □ shows the average value of the cosine of the an¬ 
gle between the orientation of the satellite and the direction of 
the centrals as a function of the separation rescaled by the av¬ 
erage virial radius (7*200). The mass of the subhaloes in the ori¬ 
entation sample (+) is kept fixed whereas the masses of the cen¬ 
tral haloes (g) are chosen to have similar or higher masses. Val¬ 
ues of (cos fi) that are significantly greater than 0.5 indicate that 
the orientation of satellites galaxies are preferentially aligned to¬ 
wards the direction of central galaxies. As for the previous subsec¬ 
tion, we use two different mass bins taken from two simulations: 
11.3 < logio(M sub /[/i _1 M©]) < 12.6 from EAGLE L100 (left) 
and 12.6 < logio(M su b/[/i _1 M©]) < 13.7 from cosmo-OWLS 
L200 (right). The line thickness is proportional to the subhalo mass 
of the position (g) sample. In both panels the physical separations 
between the pairs are normalised by the (7*200) of the haloes host¬ 
ing the central galaxies. 

The major axes of satellite galaxies, when all stars are consid¬ 
ered, are significantly aligned towards the direction of the centrals 
within their virial radius. The strength of the alignment declines 
very rapidly with radius and is very small outside the virial radius. 
There is only a weak dependence on the central subhalo mass. 

We note that by considering only stars in 7*^3f the trends 
shown in Fig.[6]and in Fig.[7]are weakened (not shown). This results 
in a less significant alignment for galaxies hosted by subhaloes with 
masses 11.3 < logio(M su b/[ h~ x M©]) < 12.6 from the EAGLE 
LI00 simulation, whereas a still significant alignment is found for 
galaxies with 12.6 < logio(M su b/[ h~ 1 M©]) < 13.7 from the 
cosmo-OWLS L200 simulation. 


4 TOWARDS OBSERVATIONS OF 

ORIENTATION-DIRECTION GALAXY ALIGNMENT 


In this subsection we report results for observationally accessible 
proxies for the orientation-direction alignment, which depend on 
the shape of galaxies as well as on their orientation, making them 
tightly connected to cosmic shear studies. All the relevant quanti¬ 
ties for the following analysis are defined in a 2D space. 

Observationally, the ellipticity is decomposed into the pro¬ 
jected tangential (e+) and transverse (e x ) components with respect 
to the projected separation vector of the galaxy pair: 


e+ = |e| cos(2T>) 
e x |§ |e| sin(2<f>) 

, , _ 1 — b/a 


( 6 ) 

(7) 

( 8 ) 


1 + b/a ’ 

where T> is the position angl^l between the projected orientation 
of the galaxy and the direction of a galaxy at projected distance r p 
and b/a is the axis ratio of the projected galaxy. 


1 In this subsection, satellite galaxies do not necessarily belong to the same 
haloes that host the paired central galaxies. 


2 The symbol <E> is used to indicate an angle between vector in 2D, whereas 
the symbol <f> (see Eq.|5j indicates an angle between vectors in 3D. 
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Figure 6. Mean value of the cosine of the angle between the orientation of the stars in the central galaxy and the direction of satellite galaxies as a func¬ 
tion of the 3D galaxy separation, rescaled by the host halo qq. The central galaxies used for the left panel are taken from the EAGLE LI00 (11.3 < 
logio(M su b/[ 6 _1 Mq]) < 12.6) while in the right panel they are taken from the cosmo-OWLS L200 simulation (12.6 < logio(M su b/[ 6 _1 Mq]) < 
13.7). In both panels the error bars represent one sigma bootstrap errors. Thicker lines indicate higher mass. The satellite distribution is aligned with the central 
galaxy out to ~ IOOt^qq. For r < 1 Or200 the alignment is substantially stronger for higher-mass satellites. 




Figure 7. As for Fig.[6]but in this case the orientation is computed for the stellar distribution in satellite galaxies and the angle is measured with respect to the 
directions of central galaxies hosted by subhaloes of different masses. The alignment between satellites and the directions of centrals decreases with distance 
but is insensitive to the mass of the host halo. 


Then the function e g+ is defined as: 


e s +( r p) — 

*7 *j\r p 


e +(i I i) 

N pairs 


(9) 


where the index i represents a galaxy in the shape sample, whereas 
the index j represents a galaxy in the position sample. The function 
e g+ (r p ) is the average value of e+ at the projected separation r p . 

Groups and clusters of galaxies, where strong tidal torques are 
expected to align satellite galaxies toward the centre of the host’s 
gravitational potential, are ideal environments to study orientation- 
direction alignment. However, the task of measuring this alignment 


has proven to be very challenging (see iKirk et alJ 120151 . and ref¬ 
erences therein). In group and cluster environments, the measured 
quantity, e g + (see Eq. n, is the mean value of the angle between 
the projected orientation of the satellite galaxy and the direction 
of the host , multiplied by the projected ellipticity of the satellite. 
Typical values of the root mean square of galaxy shape parame¬ 
ter, e *» (1 — (6/a) 2 )/(l + (6/a) 2 ), in the set of simulatio ns 
emplo yed in this study, can be found in Fig. 5 of IVelliscig et al.l 
(20151) . Those values are in broad agreement with the observed 
noise-corrected v alues (ab out 0. 5-0.6 depending on luminosity and 
galaxy type, (e.g. ljoachimi et al.ll2013h ) when all stars in subhaloes 
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Figure 8. Values of the observationally accessible proxy for orientation-direction alignment, e g + (see Eq.[9j, as a function of the projected separation r p . 
Only pairs that are separated by less than 2.5 h~ x Mpc along the projected axis are considered for this analysis. The error bars indicate one sigma bootstrap 
errors. In the left panel the orientation of the subhalo is measured using all the stellar particles whereas in the right panel only stars within the are used, 
which greatly reduces the alignment. Thicker lines indicates higher masse s. The coloured dotted lines show the average value of e g + within the virial radius 
of the central galaxy. Observational measurements from lSifon et al.l (120151) constrained the average ellipticity to be e g + = —0.0037 d= 0.0027. 


are con sidered. However, w hen only stars within if are consid¬ 
ered, Ivefiisci^etal] d2015h found typical values for e rms of w 0 .2- 
0.3, that is a factor of 2 lower than the observed value. This suggests 
that galaxy shapes computed using stars within rJlu are rounder 
than the observed shapes, potentially leading to an underestimate 
of the e g +. To quantify this effect, we would need to analyse syn¬ 
thetic galaxy images from simulations with the shape estimator al¬ 
gorithms used in weak lensing measurements. We defer such an 
investigation to future works. 


Recent observational studies of the orientation-direction 
alignment in galaxy group s and clusters repo r ted signals consis - 
tent with ze ro alignmen t ( Chisari et al.1 l20l4 ISifon et aPl2015h . 
Specifically. ISif on et al.1 d20 1 5h used a sample of ~ 14, 000 spec¬ 
troscopically confirmed galaxy members of 90 galaxy clusters with 
median mass of logio(M2oo/[ M©]) = 14.8 and median redshift 
of z t® 0.14, sel ected as pa rt of MENeaCS (Multi-Epoch Nearby 
Cluster Survey; ISand et all 120121) and CC CP (Canadian Cluster 
Comparison Project: iHoekstra et al.1 l2012h . They constrained the 
average ellipticity, within the host virial radius, to be e g + = 
(—3.7 =t 2.7) x 10 -3 or e g + = (0.4 x ± 3.1) x 1 0" 3 depending on 
the shape estimation method employed. IChisari et al.l (2014) mea¬ 
sured galaxy alignments in 3099 photometrically-selected galaxy 
groups in the redshift range between z = 0.1 and z = 0.4 of 
masses logio(M2oo/[M©]) = 13 in SDSS Stripe 82 and con¬ 
strained the alignments to similar values as ISifon et al.1 ( 20151) . 


The left panel of Fig. [8] shows the value of e g+ calculated 
for the simulations using all the stellar particles of subhaloes for 
host masses and satellite mas ses that are roughly c omp arable to the 
range of masses explored in IChisari et al.l (120141) and ISifon et all 
(12015l) . We only consider pairs separated by less than 2.5 h~ x Mpc 
along the projection axis to confine the measurement to the typ¬ 
ical extent of massive bound structures. Within the virial radii of 
groups or clusters the statistical uncertainties are large. The aver¬ 
age values of e g + for distances smaller than the host virial radii are 
« 2 — 4 x 10 -2 with errors of ^ 0.1 — 2 x 10 -2 , indicating posi¬ 


tive alignment. We repeat the same analysis using only stars within 
r h2f ( see r ight panel of Fig. BO. In this case the average value of e g + 
for distances that are smaller than the host virial radi us is consis¬ 
tent w ith ze ro, in agreement w ith the observations of IChisari et al.1 
(120141) and ISifon et all (l2015h . Using deeper observations, in or¬ 
der to probe the lower surface brightness parts of satellite galaxies, 
could represent a way to reveal the alignment that is seen in obser¬ 
vations when all stars bounded to subhaloes are considered. 

Recently, ISingh. Mandelbaum & More] ( 2015 ) measured the 
relative alignment of SDSS-III BOSS DR11 LOWZ Luminous Red 
Galaxies (LRGs) in the redshift range 0.16 < z < 0.36 o bserved 
spectroscopically in the BOSS survey (iDawson et al J2013h . As op¬ 
posed to the case of galaxy groups and clusters, these measure¬ 
ments are obtained by integrating along the line o f sigh t between 
±100 Mpc. Furthermore. ISingh. Mandelbaum & More] (l2015h re¬ 
ported the average halo masses of those galaxies, as obtained from 
galaxy-galaxy lensing ana lysis . We perform the same measure¬ 
ments as in ISingh. Mandelbaum & Morel d2015h on our simula¬ 
tions. Given the observed halo mass (logio(M{§o an [ M©]) = 
13.2) and the line of sight integration limits, we employ the cosmo- 
OWLS L200 in this analysis. 

Fig. [9] shows the values of e g +( 7 ^) from our sim ulation to¬ 
gether with the measurements from ISingh. Mandelbaum & Morel 
(12015b . Note that we have used a halo mass bin (13 < 

logio(M s c ^b/[/i _1 M©]) < 13.5) half a magnitude wide to ob¬ 
tain statistically robust measurements (iVhaioes = 1677). As for 
the case of satellite galaxies in clusters, the agreement with obser¬ 
vational results depends strongly on the subset of stars used to com¬ 
pute the galaxy orientations. When one considers all stars bound 
to subhaloes, the values obtained for e g +(r p ) are systematically 
higher than the values in observations, whereas broad agreement is 
found when using only stars inside if- 

As noted before, when only stars within if are considered, 
simulated galaxies exhibit rounder shapes than observed. There¬ 
fore, the results presented here may underestimate the values of 
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Figure 10. Diagram of the angle 'ip formed between e± of galaxy pairs at a 
distance r. 


Figure 9. Values of the observationally accessible proxy for orientation- 
direction alignment, e g + (see Eq. |9), as a function of the projected sep¬ 
aration, r p , form simulations considering all stars bound to the subhalo 
(green curves) and only stars within (mag enta curve ). The data points 
are the observational results from lSingh Mandelbaum & Morel 12015b for 
the LOWZ sample of LRG galaxies (their Fig. 19, note that in their work 
e g + is denoted (7) for its direct connection with the shear). We consider 
pairs of galaxies that have a separation along the projected axis smaller than 
100 h~ x Mpc. The error bars on the curves indicate one sigma bootstrap 
errors. When we consider all the stars, the predicted alignment is stronger 
than observed. However, when we only use stars in the part of the galaxy 
that might typically be observed, we find good agreement with the data. 


e g+ (r p ). Thus, when more observationally motivated algorithms 
would be employed to analyse the simulations, it is not guaranteed 
that the agreement found in Fig. [9] would still hold. 


5 ORIENTATION-ORIENTATION ALIGNMENT 

In this section, we present the results for the II (intrinsic-intrinsic) 
term of the intrinsic alignment that is given by the angle between 
the orientations of different haloes. We define 'ip as 

ip(\f\) — arccos(|ei(x) • e\(x + r )|). (10) 

where e± are the major eigenvectors of the 3D stellar distributions 
of a pair of galaxies separated by a 3D distance r — |r| (see 
Fig.[foJ. 

Fig. HU shows the average value of the cosine of the angle ip 
for pairs of subhaloes with similar masses at a given 3D separation 
r (in h~ x Mpc). Values are shown for four different choices of sub¬ 
halo mass, where each mass bin is taken from a different simulation 
(see legend). To estimate the errors, we bootstrap the shape sample 
100 times and take as 1-sigma error bars the 16th and the 84th per¬ 
centile of the bootstrap distribution. Values of (cos ^p) equal to 0.5 
indicate a random distribution of galaxy orientations, whereas val¬ 
ues of (cos ip) higher than 0.5 indicate that on average galaxies are 
preferentially oriented in the same direction. 

The alignment between the orientation of the stellar distri¬ 
bution decreases with distance and increases with subhalo mass. 
Comparing with Fig.|2] the orientation-orientation alignment is sys¬ 
tematically lower than the orientation-direction angle alignment. 


Beyond 50/i -1 Mpc the alignment is consistent with a random 
distribution, whereas in the orientation-direction case a positive 
alignment was found for scales up to 100 h~ x Mpc. This is sug¬ 
gestive of the direction of nearby galaxies as being the main driver 
of the orientation-orientation alignment, as a weaker orientation- 
orientation alignment naturally stems from the dilution of the 
orientation-direction alignment. 

Similarly to e g+ (in Eq. we can define the projected 
orientation-orientation e++ as: 


«++(»>) = 

I r p 


47 U 1 0 

Apairs 


( 11 ) 


where e+ is defined in Eq. [6] Galaxies are selected to have at least 
300 star particles. 

Fig.[l2lshows the projected orientation-orientation alignment, 
e++, for the same halo mass bin and integration limits as em¬ 
ployed in Fig. [9] Green and magenta curves refer to the cases where 
one uses all stellar particles in subhaloes and only stellar parti¬ 
cles confined within r^f, respectively. For comparison, e g+ (r p ) 
is overplotted in grey. As expected, the e++(r p ) profile has an 
overall lower normalization. Interestingly, e++(r p ) is steeper than 
€ s+( r P )» although the significance of this trend is diminished by 
the noisy behaviour of the e ++ (r p ) profile. 

The presence of a non-vanishing e++(r p ) profile reveals a 
net alignment of galaxies with the orientations of nearby galax¬ 
ies, thus suggesting a potential II term in cosmic shear mea¬ 
surements for galaxies residing in haloes with masses 13 < 
iQgiotMSiVlfc^M©]) < 13.5. 


6 CONCLUSIONS 


This paper reports the results of a systematic study of 
the orientation-direction and orientation-orientatio n alignment 
of g alaxies in the EAG LE (Schaveetal. 120151: ICrain et a! 


2015 ) and cosmo-OWLS dLe Brun et al.l 12014 iMcCarthv et al 


20141) hydro-cosmological simulations. The combination of these 
state-of-the-art hydro-cosmological simulations enables us to 
span four orders of magnitude in subhalo mass (10.7 ^ 

logio(M su b/[/i -1 M©]) ^ 15) and a wide range of galaxy sep¬ 
arations (— 1 < logio(r/[/i -1 Mpc]) ^ 2). For the orientation- 
direction alignment we define the galaxy orientation to be the ma- 
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r [Mpc/h] 

Figure 11. Mean value of the cosine of the angle ip between the major axes 
of the stellar distributions of subhaloes as a function of their 3D separation. 
Each mass bin is taken from a different simulation. The minimum subhalo 
mass in every bin ensures that only haloes with more than 300 stellar par¬ 
ticles are selected. Orientations are computed using all stars bound to the 
subhaloes. The orientation-orientation alignment decreases with distances 
and increases with mass. It is weaker than the orientation-direction align¬ 
ment (cf. left panel of Fig. E)- 


jor eigenvector of the inertia tensor of the distribution of stars in 
the subhalo, e±. We then compute the mean values of the angle 
cp between ei and the normalized separation vector, d, towards a 
neighbouring galaxy at the distance r, for galaxies in different sub¬ 
halo mass bins. In the case of orientation-orientation alignment, we 
compute the mean value of ip, the angle between the major axes e\ 
of galaxy pairs separated by a distance r. 

Our key findings are: 

• Subhalo mass affects the strength of the orientation-direction 
alignment of galaxies for separations up to tens of Mpc, but for 
distances greater than approximately ten times the subhalo radius 
the dependence on mass becomes insignificant. The strength of the 
signal is consistent with no orientation-direction alignment for sep¬ 
arations >> 100 times the subhalo radius (Figs. [2ll3l . 

• The difference between the orientation-direction alignment 
obtained using the dark matter, all the stars or the stars within 

to define galaxy orientations, could account for the common find¬ 
ings reported in the literature of galaxy alignment being systemati¬ 
cally stronger in simulations than reported by observational studies 
(Fig.®. Since observations are limited to the shape and orientation 
of the region of a galaxy above a limit surface brightness, simula¬ 
tions have to employ proxies for the extent of this region. 

• At a fixed mass, subhaloes hosting more aspherical or prolate 
stellar distributions show stronger orientation-direction alignment 
(Fig.0. 

• The distribution of satellites is significantly aligned with the 
orientation of the central galaxy for separations up to 100 times 
the virial radius of the host halo ( 7*200 )> within 10 7*200 higher-mass 
satellites show substantially stronger alignment (Fig. [6]). 

• Satellites are radially aligned towards the directions of the 
centrals. The strength of the alignment of satellites decreases with 
radius but is insensitive to the mass of the host halo (Fig. [7]). 



r p [Mpc/h] 

Figure 12. Dependence of (Eq. |TT}, a measure of orientation- 
orientation alignment, obtained from the simulations using an integration 
limit of 100 h~ x Mpc. Both the values for the whole stellar distribution (in 
green) and for the stars within are shown (in purple). The error bars 
indicates one sigma bootstrap errors. The results for e g + (grey curve) are 
shown for comparison. 

• Predictions for the radial profile of the projected orientation- 
direction alignment of galaxies, e g +(r p ), depend on the subset 
of stars used to measure galaxy orientations. When only stars 
within rfjf are used, we find agreement between result s from 
our simulati ons and recent obser vations from ISifon et all ( 2015| ) 
and lSingh. Mandelbaum & More] ( 20_15 )(see Figs. [8] and [9] respec¬ 
tively). 

• Predictions for the radial profile of the orientation-orientation 
alignment of galaxies, e ++ (r p ), are systematically lower than those 
for the orientation-direction alignment, e g +(r p ), and have a steeper 
radial dependence (Figs.QT]and[T2}. Although low, the non vanish¬ 
ing e++(r p ) profile reveals a net alignment of galaxies with the ori¬ 
entations of nearby galaxies, thus suggesting a potential intrinsic- 
intrinsic term in cosmic shear measurements for galaxies residing 
in haloes with masses 13 < \ogio(M^/[ h~ x M©]) < 13.5. 

For a direct comparison with the observations, in order to val¬ 
idate the models or to explain the observations, particular care has 
to be taken to compare the same quantities in simulations and ob¬ 
servations. A future development of this work will be to extend the 
comparison with observations further by using the same selection 
criteria for luminosity, colour, and morphology in the simulations 
and in the observations. 

The strength of galaxy alignments depends strongly on the 
subset of stars that are used to measure the orientations of galax¬ 
ies and it is always weaker than the alignment of the dark matter 
components. Thus, alignment models that use halo orientation as a 
direct proxy for galaxy orientation will overestimate the impact of 
intrinsic galaxy alignments on weak lensing analyses. 
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